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1. Introduction 

The investigation of the quantum mechanical behaviour of systems with classically 
chaotic counterparts is a rich and exciting field of study. Many beautiful results have 
been found over the years regarding the semiclassical quantization of nonintegrable 
systems ^ |^ and "quantum chaology" A major component of the approach 

involves the study of eigenvalue statistics, eigenfunction "scarring" etc. Classical chaotic 
behaviour, however, by definition a dynamical phenomenon, seems to be suppressed 
at the quantum level 0, §]. Apart from a brief time, during which its behaviour 
effectively mimics the classically chaotic dynamics, the quantized version of a bounded 



and conservative classical system is quasiperiodic and, therefore, predictable P, pi)|. 

Chaotic systems thus test our understanding of the relationship between classical 
and quantum mechanics. If classical behaviour is not a naturally emergent property of 
such systems in the absence of explicit measurements then it is natural to consider the 



effect of environments on the quantum system |TT|. Quantum mechanical interference 
and superpositions, the cause of the aforementioned discrepancy, can be shown to be 
destroyed on extremely short timescales, in tractable model systems, by opening these 



systems to environmental influences [12, 13|. Some features of classical behaviour can 



then emerge without having to take the semi-classical limit. 

Zurek and Paz |T^ have used the exactly solvable, inverted harmonic oscillator 



to model heuristically unstable chaotic motion near a hyperbolic point, although this 
system is certainly not chaotic itself. The Hamiltonian for this model is 

H{x,p) = P^-^. (1) 

The effect of the the environment on this system might be expected to provide clues 
to the behaviour that might be observed if we were to open a chaotic system in a 
similar manner. The oscillator can be described by a master equation in the position 
basis. Choosing the bath as a collection of harmonic oscillators in thermal equilibrium 
at high temperature, with an Ohmic spectral density proportional to a small dissipation 
parameter, simplifies this equation enormously. Analysing the time evolution of an 
initial minimum uncertainty state using this equation, Zurek and Paz have found that 
the von Neumann entropy, S{t), quickly begins to increase linearly in time with a rate 
set by the quantity, A in equation ([1|), 

S{t) ^hQ = X. (2) 



Recently, however, we have generalized this result [|T5| by considering the entropy 



production for all temperatures and for non-negligible dissipation parameters, 7^^. We 
found that equation (□) is modified to 



5(t) = /.g = JA^ + tLs-T..., (3) 
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but reduces to it when 7^4^^ ^ A, as required. 

As A is loosely analagous to the Lyapunov exponent of a chaotic system, Zurek and 
Paz conjecture that the entropy production rate might differentiate between chaotic and 
regular systems, and that this rate will also be determined by the (positive) Lyapunov 
exponent of the classical system. From the Alekseev-Brudno theorem (see p. 4 in 
also ^) 

ltoM = ^,„ = l:A^ (4) 

where I{t) is the Shannon entropy, hxs is the Kolmogorov- Sinai entropy - a measure 



of the average rate at which a classical system will lose information [|16| - and the sum 



is over all the local positive Lyapunov exponents, {A"^}, in phase space. From this 
definition, of course, it is clear that hxs is actually an entropy rate. Consequently we 
will modify the Zurek and Paz conjecture to the statement that S{t) increases linearly 
with hxs- This has the attractive feature of dealing both in the classical and quantum 
mechanical cases with the change of entropy in time. Furthermore, the equality in 
equation (^ has been replaced with the weaker statement of linear increase. This is the 
conjecture that we will consider throughout the rest of our discussion. 

The idea that quantum systems can be defined as "chaotic" or "regular" , depending 
on their behaviour under a perturbed evolution, is appealing and an old one An 
approach similar in spirit is due to Schack and Caves |]T^ , who have recently considered 
perturbed quantum maps in an information theoretic framework and have established 
a criterion for chaos known as "hypersensitivity to perturbation" . In this paper we will 
study the prototypical, theoretically important [0, and experimentally realizable 
pin , pT| quantum kicked rotor. However, we will consider the dynamics under the 
influence of a tractable, model environment and seek to find evidence for the modified 
conjecture when the parameters of the combined system are compatible with those 
assumed during the framing of the conjecture. 

The plan for the paper is as follows. In Section |^ we shall briefly outline the features 
of the classical and quantum kicked rotor relevant to the present study. In Section || 
we will derive the one-step propagators for both the (reduced) density matrix of the 
rotor and its associated Wigner function representation. The latter result allows us to 
see clearly that by opening the system we still have a quantum analogue of the classical 
kicked rotor, or rather its standard map. Also in this section we endeavour to quantify 
the effect of the environment on the coherence of the system. 

The results from our numerical analysis follow in Section ^ wherein we will produce 
firm evidence in favour of the use of entropy production as a criterion for quantum 
chaos and of the modified conjecture, at least in the parameter ranges considered. We 
will also suggest explanations for those features found not to be in accordance with the 
conjecture. The conclusions follow in Section H. 
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2. The Kicked Rotor 

The kicked rotor model has been used extensively in studies of classical chaos and its 



quantum manifestations [0, |TT|, 0, |T9|, Consequently the literature is vast, but 



for convenience we outline here those features of its complex behaviour necessary for 
an understanding of our model and results. Even for readers familiar with the model a 
perusal would be useful in order to establish the notation used. 

2.1. Classical Features 

The system Hamiltonian with which we will deal is that of a periodically kicked rotor, 
HrotoT = y ~ ^ cos(27rg) ^ 5{t-n), (5) 



where p is an angular momentum and g is a scaled angle variable. Integrating the 
equations of motion from immediately after a delta function "kick" to just after the 
next such kick will give a version of the famous standard map 

K 

Pn+i =Pn- TT 27rg„+i, 

2tt (6) 

qn+l = iqn+Pn) (mod 1). 

The map is periodic, with period 1, in both directions of phase space - in the angle 
direction by definition and in the momentum direction due to the form of the map. 
This means that one can simply view trajectories projected back onto the unit square, 
[0, 1) X [0, 1), in order to study the dynamics for different initial conditions. However, 
even though there is periodicity in the momentum direction, the momentum is actually 
defined on the infinite interval [— cxd,cxd]. The phase space has, therefore, the topology 
of a cylinder. 

The nonlinearity of the map is controlled by strength of the kick, K. The trivial, 
regular motion when K = 0, becomes more and more complicated as K is increased 



||T8| , IS]. When K < Kc ^ 0.97... some KAM tori span the cylinder and prevent 
diffusion in the momentum direction. As K increases still further these tori are 
destroyed, allowing the momentum to diffuse, i.e. 

Mn) ~ p{0)f) = D{K)n, (7) 

at a rate determined by the diffusion constant, D{K). Here, the averaging is over a 
distribution of initial conditions. The diffusion constant can be calculated analytically 



23], and it is found that correlations in the angle variable from kick to kick 
influence its form considerably. If we ignore such correlations - an increasingly accurate 
approximation as K — oo - the angle variable can be assumed to be randomly but 
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uniformly distributed on the unit interval. This quasilinear approximation allows us to 
write 

0,.iK)^'^{^)\ (8) 
Inclusion of these dynamical correlations leads to corrections which decrease in 



magnitude as K gets larger, i.e. it is found |T9|, p3 



where J2 denotes a Bessel function of order 2. 

For all finite values of K, KAM tori remain around elliptic points in phase space, 
of which there are an infinite number. These stable islands, surrounded by a chaotic 
sea, characterize what is known as a "mixed" phase space, meaning that two nearby 
trajectories either remain "near" to one another for all time or diverge exponentially, 
depending on whether they are to be found in a regular region or a chaotic one, 
respectively (See figures and [§]). Hence the Lyapunov exponents for the system 
are dependent on the starting point in phase space, but gradually converge to a uniform 
value as ^ cxD. 



From Pesin's theorem |24|, the kicked rotor, being a chaotic system with one degree- 



of-freedom, possesses a positive KS entropy, hxs, which will be equal to its positive 



Lyapunov exponent for K sufficiently large. Chirikov has calculated this quantity 
and found 

hKS ^ In (y) . (10) 
2.2. Quantum Features 

The periodicity of the angle variable gives rise to a numerable and discrete set of 
momentum eigenvalues defined by 

p\l)=pi\l)=27rhl\l), (11) 

where / is an integer. The corresponding eigenstates are obviously eigenstates too of the 
kinetic energy operator Hq = p'^/2, and we write 

Hq\1) = Ei\l) = 2n^hH^\l). (12) 

The Hamiltonian of equation (^ leads to a one-step unitary transformation 

U{tn+l,tr,) = UkUf, (13) 

which defines the quantum standard map. Here, tn is the time after the nth kick and Uk 
is the unitary kick evolution operator. It can be written in the momentum representation 

as 

{l'\Uk\l)=^''-'Jl'-l{K/4n^h), (14) 
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with Jm denoting a Bessel function of integer order m. Uf is the operator corresponding 
to unitary free rotation between the deha function kicks, and has the momentum 
representation 

{l'\Uf\l) = exp{-i27T^hP) 5n. (15) 

The unitary dynamics was first investigated by Casati at al. |jlO[ using initial states 
highly localized in momentum. It was found that {{p{n) — p{0))^) shows diffusive 
behaviour in approximate agreement with the classical case, but only up to a finite 
break time, n*, after which the momentum ceases to grow further and quasiperiodic 
motion around some finite average is found. This defines a localization length via 

mn)-p{0)y) = i2nhfL', (16) 

for n > n*. Good analytic and numerical evidence now exists pointing to the fact that 
the eigenstates of the unitary evolution operator of equation ( |T3| ) above are localized in 
momentum and have essentially discrete spectra on the unit circle of the complex plane. 



These assumptions lead to the following estimates for L and n* ^0 



(17) 



and 

n* ^ 2L. (18) 

It is important to note that both quantities are inversely proportional to the square of the 
scaled unit of action, h. We also note that classical correlations, which can both increase 
and decrease the diffusion constant from its quasilinear value, Dgi, (see equations (§) 
and (^) can also directly affect the quantum dynamics via this localization mechanism. 

The discrepancy between the classical and quantum diffusive behaviour of the 
kicked rotor is therefore known to be a consequence of quantum interference. If 
the phase coherence could somehow be destroyed one might expect to recover the 



classical diffusion. Various authors have considered this problem [0, |2^, |2^, and 
have examined the effect of an environment on localization. However, previously the 
relationship between the classical and quantum entropy production, the former being a 
quantitative measure of chaos, has not been studied. This approach is important since 
we have growing evidence that it is a generic one. 



3. Open Dynamics 



We now proceed to derive the propagator for the density matrix of the rotor after tracing 
out the environmental degrees of freedom, i.e. the propagator for the reduced density 
matrix. This will enable us to calculate the von Neumann entropy at each time step, 
upon iteration of an arbitrary initial state. 



3.1. The Reduced Density Matrix Propagator 

We will restrict attention to the situation between kicks initially and consider the free 
rotor to be coupled through its momentum to a thermal bath of harmonic oscillators. 
The Hamiltonian for this system, 

^ - y + ^2^, + + ^ \ ''"^'^ ^^^^ 

can be written as 

H = j + Y: ^Mbibk + h + pE^(4 + bk), (20) 

k k 



where Ck = \ Ti/ {2mkUJk)ck is a renormalized coupling and 



rUkUk f ^ , , Pk \ 



-T^T- I qk + %^ I 



qk ' 



2h \ rnkUk, 

As [qj,Pk] = if>,6jk, these annihilation and creation operators satisfy = Sjk- 

We choose an uncorrected initial product state 

p(o)=p^(o)®npL (21) 

k 

where p^f^ is the density matrix of the kth harmonic oscillator in thermal equilibrium at 
temperature T, defined by 



Pth 



1 — exp 



-1 



—nhuj 



To specify the problem further we also need to choose a physically reasonable 
environmental spectral density p9[, /(t^), defined formally by 

/M=i:.(.-..)^. (23) 

More generally, if we assume every environmental oscillator to have equal mass, we can 
write H 

JM^^M^, (24) 

with a smooth density S. In order to avoid the influence of unphysical, extremely high 
frequencies we assume the existence of an upper, or cutoff frequency, ujc, and choose a 
continuum of oscillators distributed according to the density 

S{uj)c^{uj) = 2mr7w^exp(-— ), (25) 
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with q > 1, which reduces to that used by Caldeira and Leggett [|29[ when we choose 
a so-called Ohmic environment with q = 2. Hamiltonians and initial conditions of this 
Ohmic type have been considered by Jiushu et al. in |^3[ and the exact time evolution 
for the reduced density matrix has been obtained. We proceed along similar lines but 
consider the general case with g > 1 being the only assumption. We arrive finally at an 
expression for the reduced density matrix of the rotor in the momentum representation 

{m\p'{t)\n) = exp (^~{E„. - En){t + A{q,t)) - ^^"^ ~ t)) (m|p^(0)|n). (26) 

The environmental terms in this expression, given by A{q,t) and B{q,t), have been 
evaluated explicitly for all g > 1 [0. The expressions are particularly simple, however, 
for the Ohmic case, and we concentrate on this important case from here on. We find 

A{2, t) = 2r]{aTctan{uJct) - Uct), (27) 

and 

B(2.t) = |,„(i +,i„n(i + (j^n (28) 

where the usual notation of /5 = l/^ksT) has also been used. 

There are a few points to note here. Firstly, were K = in equation (^, then 
equation (|^) would describe fully the time evolution of any initial state, p'^(O), in the 
momentum representation, for a classically nonchaotic rotor. Secondly, we return to the 
correct unitary dynamics , corresponding to conservative evolution of the rotor, when 
rj is set to 0. (This amounts to setting all coupling constants, c^, to zero.) Finally, we 
will consider equation ( ^61) to describe the time evolution of the quantum kicked rotor 
between kicks and so will set t = 1 in equations (p6D, (p7|) and (p8|). The heat bath is 
assumed to be sufficiently large to remain in thermal equilibrium as it interacts with 
the system between each kick. For this assumption to be valid the coupling constants, 
{cfc}, are taken to be small. Hence, for the spectral density in equation (|25|), r/ <^ 1, a 



common assumption 14, 30 



The system evolution, from just after the kick at t = n until just before the kick at 
t = n + 1, where n is an integer, is given by a master equation of the form 

{l'\p~{n + l)\m') =Y,G^{l',m',l,m){l\p+{n)\m), (29) 

l,Tn 

where 

(/', m, /, m) = G{.(/', m', /, m)T{l, m), (30) 
Gq{1' ,m' ,l,m) is the unitary summation kernel given by 

Gij{l', m', I, m) = exp{--{Ei - Em)) ^vi^m'm (31) 
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and 

r(/, m) = eM-kEi - E^)A{2, 1) - ^ILZP^b{2, 1)) (32) 
n n 

describes the non-unitary influence of tlie environment on tlie dynamics. 

It is assumed tliat during eacli kick botli free evolution and the influence of the 
environment are negligible. We are then led formally to the summation kernel for the 
kick, defined by the master equation 

{l'\p+{n + l)\m') = Y,G''{l',m',l,m){l\p-{n+l)\m). (33) 

Choosing now the shorthand notation used by Dittrich and Graham |3^ of k := K/An^h 
and bm{x) := i"^Jm{x) we can write explicitly 

G\l',m',l,m)=b,^i{k)bl,_^{k). (34) 



Even as the density matrix evolves, according to master equations (^Oj) and (p^) 
above, it must, of course, remain normalized, hermitian and positive These 
properties in turn impose conditions on the kernels of equations ( PO]) and ( |5^ which 
can easily be shown to hold ||3^, implying that the evolution of the reduced quantum 
statistical density operator is a valid one. 

The one time-step evolution of the density matrix is now simply the convolution 
of and G'' defined in equations (|30D and (0) above, respectively, and it is easily 
confirmed that 



{l'\p+{n + l)\m') = Y,G{l',m',l,m){l\p+{n)\m), (35) 

l,m, 

where 

G{1', m', I, m) = Gc{l', m', I, m)T{l, m), (36) 

in which 

Gc{l',m',l,m) = {l'\U\l) {m\U^m') (37) 

is the unperturbed, unitary propagator. When there is no environment, i.e. as r/ — > 0, 
then r(/, m) — > 1 for each / and m. 



3.2. Destruction of Coherence 

It was stated above that the term r(/,m) in equation ( p6D makes a non-unitary 
contribution to the dynamics. In this section we would like to quantify this statement. 

A measure of decoherence is the decay probability. Pi, of a quasienergy eigenstate in 
a single timestep. In order to calculate Pi we first consider the probability of transition. 
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P^{s\r), from an quasienergy eigenstate |r) to a quasienergy eigenstate |s) in a single 
timestep. In fact, if p^(0) = \r){r\ then we may write 

pi(s|r) = (s|p+(l)|s). (38) 

The propagator in the quasienergy representation is therefore required, and the master 
equation ( P5| ) maybe be rewritten 

(ai|p+(n + l)|/5i) = J2 G'«'^(«i,/?i,«o,/5o)(«o|p+(n)|/5o), (39) 

ao,/3o 

where 

U\ao) = e''^-o\a,), (40) 
and similarly for and One easily finds 

G"'^(ai,/3i,ao,/3o) = Ee*^'^"^"'^'^^(«ilOa|ao)(/?o|m)(m|A)r(/,m). (41) 

l,rn 

Thus, because of equations (^) and (|39D, we find 
P\s\r) = G^"(s,s,r,r) 

= E(s|/)(/|r)(r|m)(m|s)r(/,m). (42) 

l,m 

We can now use this expression to calculate Pi by summing over all possible final 
states, \s) (different from |r)) and averaging over the various possible initial states (see. 



for example, reference |2 



A = ^ E ^'(^1^). (43) 



D ^ 

where D is the dimension of the truncated basis (See Section ^). Inserting equation (^) 
and using the fact that 

J2 {s\l){m\s) = 6^,1 - {r\l){m\r) (44) 

we find 

^i = l--^E7., (45) 

where we have defined 

7. = EK^I0ri(r|m)rr(/,m). (46) 

l,m 

We may estimate this expression analytically if we assume that every 

quasienergy eigenfunction is exponentially localized according to 

|(r|OP^expf-^^V (47) 
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Here we have denoted the centre of localization by and L is the localization length 



defined previously in equation (|T7|). However, not every such eigenstate has this 
structure, as certain "double-hump" quasienergy eigenstates exist, so-called because 
they have two localization centres. 

For accuracy therefore, we have calculated Pi numerically and examined the effect 
changing K, % and /?, the environmental inverse temperature. We present our results 
in Section |. 

3. 3. The Wigner Distribution Propagator 

We would like now to examine the classical limit of our system. To do so we choose to 
work in the Wigner representation for a phase space with the topology of a cylinder. 
The notation in this section is exactly that used in refs. [^, ^ , in the latter of which 
this representation is succinctly descibed. 

The Wigner distribution on our cylindrical phase space at time n, Wn{q,p), consists 
of 5-functions concentrated on the definite values of momentum, p = ttM, for / an 
integer. One therefore writes 

Wniq,p) = E Wi'\q)Sip - Tihl) (48) 
I 

for the Wigner distribution at this time. Its (one) time-step generator, Gw, is defined 
by 

wilM = E ['dqGw{lq,l,q)WP{q), (49) 

and can be expressed in terms of the one-step density matrix propagator, most 
conveniently in the momentum representation, using 

GwiJ, q,l,q)= E E ^i,rm+m2^i,h+i2 exp{27ri[(mi - ma)? - (h - /s)?]} 

xG'(mi,m2,/i,/2)- (50) 

It will be convenient to make the transformation from the density matrix propagator 
to the Wigner function propagator for both the free and the kicked part separately. A 
convolution of the two will then yield the desired quantized map in the Wigner function 
representation. 

Use of equation (|30D in equation (|50|) gives 

G^^.(I,g,Z,g) = 5|,Ga(x,A,l), (51) 



0(J 



fc=— oo 

i.e., proportional to a periodic Gaussian with period 1. Here, x and A are defined by 
X = q~ q-Tini{l + A{1)), (52) 
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A = 2hB{l). (53) 

This periodic Gaussian reduces to a periodic (5-function with the argument q — q — ttM 
for ?7 — > 0, as expected. Similarly, use of equation (|3^) in equation (|^) gives 



Gwil q, I q) = i'"' exp{27ril(g - q)}J-_pk sin(27rg)) ^ exp{-47rim(g -q)}. (54) 

m 

There remains only to convolute the two propagators to give 

Gw{l,^,l,q) =J2 dqG'^{l,^,l,q)G{y{l,q,l,q) 
1 ^° 

= J=_^(2A;sin(27rg)) Ga(g - g - 7rhl{l + A{1)), 2hB{l), 1). (55) 



When ?7 -C 1 and h this quantum map can be shown to reduce to the 
classical map of equation (|]), as required. 

4. Results 

We now present the results of our numerical calculations. For each run we began with 
an initial coherent state centred at some point {p, q) in phase space. We choose Planck's 
constant, /i to be a rational multiple of 

7 (56) 

This avoids any resonances which might occur between the frequency of the 

driving force and those frequencies associated with the energy levels, {-En}, of the 
unperturbed, regular rotor for rational h. (Such resonances give rise to a superdiffusive, 
quadratic increase of the rotor's energy with time in contrast to the linear increase 
found for the classical rotor.) The infinite dimensionality of the Hilbert space forces 
us to truncate the momentum basis in which we work. We have chosen, for the 
generation of the results presented below, to use a truncated basis set which is 
symmetric about the zero momentum eigenstate. That is, we have used the set 
{| — n), I — n + 1), . . . , I — 1), |0), I + 1), . . . , |n — 1), |^)}, where n was chosen large 
enough to keep the density matrix normalized to an accuracy of better than 1 part in 
10^; i.e. 

Tr[p5(t)] = l, (57) 



for the time of interest. Typically, n ranged from ^ 250 up to ~ 450, depending on the 
parameter values used. 
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A propagation through one time step was then effected by computing equation (p^) 
followed by equation (|33|). At each stage the von Neumann entropy, S(t) was calculated 
by using the approximation 

n 

Sit) := -TT[ps{t)\npsit)] ~ - E ^^l^A^, (58) 

i=—n 

where {Ai| — n < i < n} is the set of eigenvalues of the density matrix ps{t) at that 
time. The entropy is a maximum if and when each Aj = l/{2n + 1), giving an upper 
limit to the entropy of Smax = ln(2n + 1) as a result of the truncation. 

Finally, as we wish to investigate the modified conjecture it is reasonable to make 
the same characterisation for our environment as Zurek and Paz |T^. The spectral 
density is Ohmic but a high temperature environment is also required. Consequently, 
the parameters are chosen so that ksT > huc- 

4-1. Qualitative Comparisons 

4-1.1. Local KS Entropy When the nonlinearity parameter, K., of our rotor is greater 
than or equal to 6 there are no discernable stable islands in phase space plots of the 
classical standard map. They do still exist, of course, but the fraction of the unit square 
they occupy at these values is small. In such a situation the variation in the Lyapunov 
exponent will be small throughout phase space and the inverted oscillator model will 
be more likely to be considered adequate. However, when the value of K is so small 
as to allow KAM islands to exist and occupy a significant fraction of phase space we 
must expect the local Lyapunov exponents to vary considerably. All of these local 
Lyapunov exponents contribute to the sum in the calculation of the global KS entropy 
in equation @). However, we can also define a local KS entropy by restricting the 
sum over local Lyapunov exponents to those corresponding to the points in an initial 
classical probability distribution. Thus, in view of the aim of this investigation, we 
might expect that the rate of von Neumann entropy increase will depend strongly upon 
precisely where in a mixed space we locate an initial coherent state, i.e. upon the local 
KS entropy corresponding to its classically analogous initial probability distribution. 

To illustrate the usefulness of the local KS entropy concept we considered first of 
all the time evolution of two initially Gaussian, symmetric probability distributions 
in phase space, one centred in a classically regular region when K = 1, with the other 
centred in a chaotic region for the same K value. The variance of these distributions was 
chosen as h/2 in order to render them classical analogues of the minimum uncertainty 
quantum initial states to which their evolutions will be compared. 

In figures and [0] we plot the result of taking 100 sample points from each of these 
distributions and plotting the first 300 points of the classical trajectories corresponding 
to each such point. The regularity of the dynamics in figure where the initial state is 
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centred at (0.0, Oh), h being Planck's constant as before, can be compared to the chaotic 
trajectories in figure resulting from the initial state being centred at (0.5, 2/i). 

Our next step is to calculate the KS entropy corresponding to each of the points 
used to sample the Gaussians and then to take the average. This we define as the local 
KS entropy corresponding to the initial states. The calculations are done for a range of 
K values from to 6 giving rise to progressively less mixed phase spaces. In figure P 
we see that this local KS entropy differs greatly for these two states when 1 < < 4, 
i.e. after a significant number of invariant KAM curves have been destroyed and before 
the K value at which the KAM island surrounding the elliptic fixed point at the origin 
disappears completely from phase space plots |]18|, As K increases further these 
local KS entropies agree more and more, and can be well approximated by the global 
KS entropy 12^ of equation (p!oD. 



4- 1.2. Local Quantum-Classical Entropy Correspondence Though global KS entropy 
can, of course, be defined for all values of K, it must be numerically determined when a 
largely mixed phase space prevails. We show here, however, that it is not the quantity 
to be used in any prediction of the quantum mechanical entropy production rate. The 
correct quantity is the local KS entropy. 

To show this we next considered the quantum mechanical time evolution of the two 
coherent states analagous to the two classical distibutions used above. For each state 
we have calculated the von Neumann entropy production rate over 20 time steps and 
for the set K = 1, 2, . . . , 5, values at which the locally calculated KS entropies differ 
greatly (See figure 0]). The results corresponding to the initial, classically regular state 
centred at (0.0, Oh) can be seen in figure and those corresponding to the (classically 
chaotic) state centred at (0.5, 2h) can be seen in figure p. In the former case the entropy 
production rates refiect very well the local KS entropy values of figure . The initially 
low values of each, when K = 1,2, increase when K = 3, once more when K = 4, and 
finally jump dramatically when K = 5. However, the situation for the coherent state 
situated initially in a classically chaotic region of phase space is strikingly different. 
Indeed, the steady increase in plateau height with K mirrors well the local KS entropy 
of figure as it, too, increases steadily with K. No sudden jumps are apparent, and 
only when = 5 do the entropy production rates of both initial states approach one 
another, just as they do classically. 



4.1.3. Overlap Effects in a Mixed Phase Space Entropy production, thus, does depend 
on where we initially place our coherent state but we expect this variability to decrease 
with increasing values of K as the classical phase space becomes more uniformly 
chaotic. However, the initial states for = 1 we have considered above have been 
well localized (figures and 0) in either entirely regular or entirely chaotic regions. 
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with no significant overlap with nearby regions which may have very different classical 
behaviour. In order to investigate the effect significant overlaps have on the rate of 
entropy production we set K = 4, where, as can be seen from the classical phase space 
plot of figure 0, one large regular region still exists in a sea of chaotic trajectories. If we 
locate our initial coherent states along the diagonal which intercepts this KAM island, 
with their centres equally spaced at (0,0), {h, h), {2h,2h), . . ., we find an interesting 
feature. The rate of entropy increase is lower when the initial state overlaps the regular 
region and only when a significant overlap exists initially with the chaotic sea does the 
entropy begin to increase at the robust rate which is found for all the states with little 
or no overlap with this island. Classically we see analagous behaviour once more if we 
examine the locally averaged KS entropy for initial distributions which correspond in 
their mean and variance to those quantum states of figure [0]. We show the results of 
this calculation in figure Here, as in the quantum case, an asymptotic approach to 
a constant KS entropy value (of ~ 0.8) is found as the initial distibution overlaps to a 
lesser degree with the stable KAM island. 

It seems, therefore, that entropy production is indeed a suitable criterion by which 
to distinguish between regular and chaotic motion. Recently Zarum and Sarkar 
produced qualitative evidence in support of this claim for the kicked rotor model. They 
quantized on a torus and perturbed the rotor in a different, but noisy and non- unitary 
way, leading to an entropy increase. Initial coherent states uniformly distributed on the 
torus gave rise to a corresponding distribution of (constant) entropy production rates 
which reflect in a remarkable way the classical KS entropy profile for initial points in a 
mixed phase space. Similar results have recently been obtained for the quantum kicked 
top model [ P7[ . 

4-2. Quantitative Comparisons 

The quantitative analysis of the relation between hq and hxs has been the goal of this 
paper and we now turn our attention to it. 

4-2.1. Temperature-induced Correspondence Firstly we examine the effect that 
increasing the environmental temperature has on quantum-classical correspondence. To 
do so we have once more iterated in time a normally distributed initial classical state 
consisting of 25000 points. The distribution was chosen to be centred at the same point, 
and to have the same variance, as a minimum uncertainty quantum state centred at 
(0.5, h). We have chosen /i = O.I7 and i^' = 3 to examine the correspondence in time 
between the classical average value and the quantum expectation value of the squared 
momentum variable in figure ^j. Clearly, the classical behaviour, marked with stars 
(*), is approached to a greater degree as the temperature is increased (equivalently, as 
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(3 = l/ksT is decreased) despite the fact that h remains fixed. 

Now, in figure |1^ we examine, for the same parameters, the consequences for 



the rate of von Neumann entropy production and in figure the consequences 
for the von Neumann entropy itself. The truncated basis chosen for these runs as 
{\ — n), \ — n + 1), . . . , \ — 1), |0), I + 1), . . . , |n — 1), |n)}, where n = 250. This imphes a 
maximum entropy value of In 501 ~ 6.2. 

For inverse temperatures (3 = 1.00, 0.10 and 0.01, and for the time regime shown, we 
find that the von Neumann entropy is not near enough to the saturation value for the 
linearity of its increase with time to be affected. Eventually, however, the entropy will 
begin to saturate even for these low temperatures and the rates of entropy increase will 
fall below that of their respective "plateaux" shown in figure [|l^]. A feature noi expected 
on the basis of the inverted oscillator model [l^, jTS), and yet clearly in evidence here, 
is the dependence on temperature of the plateau height. Hence the need for a modified 
conjecture to take this effect into account. 

When we increase the temperature still further, so that (3 = 0.001, we note 
two generic features. Firstly, the time at which the approximately constant entropy 
production is attained, i.e. the time at which the plateau begins, decreases with the 
increase in temperature. This is due to the greater effectiveness of high-temperature 
environments in destroying quantum coherence. Secondly, we note that the plateau 
exists for a shorter length of time as the temperature is increased. This is because the 
increase in the plateau height due to the increase in temperature causes the entropy 
to approach its saturation value at an earlier time (by t = 8 here). We note that the 
first of these features was mentioned by Zurek and Paz in the framing of their original 
conjecture They also expected a finite lifetime for the plateaux due to saturation 
effects. However, the increase of plateau height with temperature seen here shortens 
this lifetime still further, making the plateaux seem more and more like maxima as the 
temperature is increased. 

For the highest temperature chosen, with (3 = 0.0001, we see that the initial 
coherence is lost so quickly as to render the dynamical localization mechanism ineffective. 
Indeed, a linear increase of (p2(t))/2 with the time, t, is seen in figure at this 
temperature, albeit at a lower rate than the classical value. This discrepancy, however, 
disappears as the temperature is increased still further. With regard to entropy and its 
production rate we see the above-mentioned features strikingly displayed: the entropy 
production rate has a maximum rather than a plateau; this maximum is greater than 
the plateaux heights of lower temperatures, and the maximum is reached at the shorter 
time of t = 3. 



4-2.2. Factors Affecting the Decay Probability Evidence that increasing one or more 
of the temperature, h and K, the magnitude of the kick, destroys the quantum coherence 
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more efficiently can be seen in figures [|T2| and [|1^]. There we have plotted the single-step 
decay probability, Pi, described in section |^ against for different temperatures (figure 
||T^ ) and values of the quantum of action (figure |13[)- These figures clearly indicate 
that to destroy the dynamical localization mechanism more efficiently one should change 
these parameters in this way in order to reduce the coherence time, ric. Furthermore, 
by increasing K and decreasing h we also increase the value of the break time, n*, 
according to equation (|T^). This gives the environment "more time" to act as the 
situation n* ^ rif. develops. 

These observations are important in any consideration of the rate at which the 
system produces entropy. Entropy production in the system is a manifestation of the 
increase in phase space area (or volume for higher dimensional systems) that a wave 
packet occupies when it is perturbed by an environment. Localization, therefore, and 
the consequent restrictions this places on the extent to which this area can increase, 
must be considered as a factor which can affect the entropy. We note that the inverted 
oscillator model, notwithstanding its usefulness as a guide to intuition, cannot take this 
effect into account. 



4-2.3. The Modified Conjecture In order to examine the validity or otherwise of the 
modified conjecture we plot, in figures to |T^, the constant rate of von Neumann 
entropy production, hg, versus the corresponding, classical, local KS entropy for given 
initial states, with different environmental parameters and values of the quantum of 
action. 

In figure [0] we see a linear relationship but with oscillations superimposed. Here, 
h = 0.057 cind f3 = 0.10. Evidence that these oscillations have their origin in the 
dynamical localization mechanism can be gathered from the fact that the values of 
K at which the plot departs most from linear behaviour are precisely those values for 
which the quasilinear approximation to the classical diffusion coefficient, D{K), is worst 
(see equation (H) and references [|19|, |2^). Thus, on account of the link between both 
the localization length, L, and the break time, n*, given by equations (|l^ and ([T8|) , 
respectively, and the classical diffusion coefficient, we see that these values of K give 
rise to exceptionally small or large localization lengths and break times. This, for the 
reason given above, affects the rate of entropy production. 

To investigate this effect further we have increased the temperature tenfold in figure 
|[T5| , i.e. f3 = 0.01. We have also chosen to start initial minimum uncertainty states at 
five different points in phase space. Since the quasienergy eigenstates form a complete 
set, each initial state can be expanded in this basis and will have different expansion 
coefficients. Differing levels of importance will then be attached to the various basis 
states for each initial condition. Since the localization length varies from eigenstate to 
eigenstate we might therefore expect to see slightly different entropy production rates 
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for each such state, but with the hnear behaviour still clearly in evidence. This is indeed 
the case. We also note that the slope of the graph has increased with the increase in 
temperature. 



Running at higher temperatures when h = O.O57 as in figures and |T5| quickly 
leads to undesirable saturation effects so we henceforth choose h = O.l'-f and examine the 
interplay between dynamical localization and the modified conjecture with this larger 
value for the quantum of action. 

In figures [|r^, [T^ and [|r^ we have done just that, but for values of the inverse 
temperature of (3 = 0.1, 0.001 and 0.0001, respectively. The starting point was chosen as 
(0.5, h). Thus, the same parameter values are used here for values of K from 3.5 to ~ 25 
as were used in figures @, flOl and ||lll for if = 3. Figure f^, especially, should give 



one an indication of the increased effectiveness of a larger environmental temperature 
in destroying the dynamical localization. This is manifest in these figures. As the 
temperature is increased we see a transition from quite a poor linear relationship in 
figure with large fluctuations due to the existence of a localization effect (especially 
for smaller K values) to obviously much better linear relationships in figure ll^ and 
The fluctuations, though still present, are much reduced at these temperatures 



since there is now little or no localization to inhibit the spread of the evolving wave 
packet in the momentum direction. Once again we see that increasing the temperature 
has the effect of raising the slope. 



5. Conclusion 

We have shown that the rate of von Neumann entropy production for an initial state 
of the open quantum kicked rotor increases approximately linearly with the local KS 
entropy of a corresponding classical distribution. This approximation becomes more 
exact as the strength of the environment becomes great enough to destroy on a short 
enough timescale those quantum coherences which lead to dynamical localization in 
the absence of an environment. However, even in that "transition" region when 
the localization is just noticably destroyed, linear behaviour is observed, albeit with 
localization-induced fluctuations superimposed. 

We have also shown using this model, but believe it to be generally true given its 
prototypical nature, that: it is generally the case that the more classically chaotic 
a system is, the greater the rate its quantum counterpart will lose information to 
a perturbing environment; b) all system and environmental parameters being equal, 
different entropy production rates can still be seen depending on where one places the 
initial state in phase space. These locations correlate excellently with the classical mixed 
phase space structure, and finally; c) the constant rate of entropy increase depends very 
much on the nature and strength of the perturbing environment. 
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We note that observation b ) here is not unique to the approach of couphng quantum 
chaotic systems to a thermal environment, as has been done here, or of otherwise 



perturbing it 0. Indeed, recently Mirbach and Korsch have introduced a 
phase space entropy-like quantity, S{p, q), which measures the (time averaged) spreading 
of mimimum uncertainty wave packets centred initially at {p,q). No environment is 
introduced in their approach. Contour plots of this quantity show increasing similarities 
to both the classical phase space structure of Poincare sections of a harmonically driven 
rotor, and to a suitably defined c/asszca/ phase space entropy, as h is decreased. However, 
no connection is made with standard quantities used in the study of chaos, such as the 
KS entropy considered here. 

These results, we speculate, provide good criteria by which to distinguish between 
the quantum mechanical counterparts of regular and chaotic classical motion in general, 
even in multi-dimensional systems. 
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Figure 1. The squares arc 100 sample points from an initial Gaussian distribution 
centred at the origin (0.0, Oh) with a variance of fi/2 corresponding to a quantum initial 
state of minimum uncertainty. We choose h = O.O57/ {2tt). The dots mark the first 300 
iterations of each such initial point, thereby giving an indication of the predominantly 
regular resultant trajectories. 
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Figure 2. The squares are 100 sample points from an initial Gaussian distribution 
centred at (0.5, 2/i) with a variance of fijl corresponding to a quantum initial state 
of minimum uncertainty. We choose fi = 0.057/(27r). The dots mark the first 300 
iterations of each such initial point, thereby giving an indication of the predominantly 
chaotic resultant trajectories. 
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Figure 3. Tlie locally averaged KS entropy of each distribution in figures ||T[ and ||^, 
corresponding to the filled squares and the filled circles, respectively, plotted against 
K. 




Figure 4. The change in the entropy production rate, hq, with time, t, for various 
values of K. The initial state corresponds to that plotted in figure giving rise to 
regular trajectories. Parameter values used are h — 0.1j,P = 0.1. 




Figure 5. The change in the entropy production rate, hg, with time, t, for various 
values of K, and using the same parameter values of figure The initial state 
corresponds to that plotted in figure 0, giving rise to chaotic trajectories. 




Figure 6. A classical phase space plot for K = 4 with the stable island clearly visible. 
The superimposed symbols correspond to those in figure Q and mark the centre of 
each initial coherent state for which the entropy rate is plotted in that figure. 




Figure 7. The rate of von Neumann entropy production for various starting points 
in the mixed phase space of figure and denoted by the same symbols here. Here, 
h = 0.057. 



1 2 3 4 5 

qO (=pO) / h 



Figure 8. The classically calculated KS entropy for various initial Gaussian 
distributions centred along the diagonal qO — pO in the mixed phase space of figure |^ . 
The position of the filled circles along the horizontal axis indicates the centre of the 
initial distribution. Again we choose h = O.O57. 
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Figure 9. Here we see that increasing the temperature of the environment 

(cquivalcntly, decreasing /?) has the effect of restoring the correspondence between 
quantum and classical expectation values. See the text for details. 




Figure 10. Here we see the consequences that restoration of the quantum-classical 
correspondence has for the von Neumann entropy production rate. See the text for 
details. 
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Figure 11. Here we see the consequences that restoration of the quantum-classical 
correspondence has for the von Neumann entropy. See the text for details. 




Figure 12. In this figure we plot the variation of the single-step decay probability, 
Pi, with K, for various different temperatures. We have chosen h = O.I7. See section 
H and section || for details. 
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Figure 13. In this figure we plot the variation of the single-step decay probability, 
Pi, with K, for two differnt values of h. We have chosen /3 = 0.01. See section |^ and 
section pi for details. 
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Figure 14. The constant entropy production rate plotted against tlie corresponding 
local KS entropy. The starting point in space is (0.5, 2/i,), where h = O.O57. Here 
(3 = 0.1. 
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Figure 15. The constant entropy production rate plotted against tiie corresponding 
local KS entropy for five different starting points in phase space. Here, as in figure 
we choose h = O.O57 but f3 = 0.01. 




Figure 16. The constant entropy production rate plotted against the corresponding 
local KS entropy. The starting point in phase space is (0.5, /i), where h = O.I7. Here 
(3 = 0.1. 
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Figure 17. The constant entropy production rate plotted against the corresponding 
local KS entropy. The starting point in phase space is (0.5, /i), where h = O.I7. Here 
f3 = 0.001. 
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Figure 18. The constant entropy production rate plotted against the corresponding 
local KS entropy. The starting point in phase space is (0.5, /i), where h = O.I7. Here 
f3 = 0.0001. 



